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Abstract 

Energy level spacing statistics are discussed for the octagonal tiling, a two 
dimensional quasiperiodic structure. A recursion relation is written for the 
probability distributions of variables defined on finite size approximants to 
this quasiperiodic tiling, using their property of similarity under inflation. 
Three types of distribution functions are introduced and determined by a 
combination of numerical and analytical techniques - these are likely to be of 
general utility in systems lacking translational invariance but with inflation 
symmetry. 
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Electrons in quasi-periodic solids are expected to have electronic properties that are in 
some sense intermediate between those of perfect crystals and disordered solids. This comes 
from noting that on the one hand, the lack of periodicity of the tilings precludes construction 
of extended Bloch states as in the case of crystals, thus reducing the possibility of metallic 
behavior or these systems. On the other hand, in a strongly disordered solid, electronic states 
existing for a given energy in a given neighborhood of the sample are located exponentially 
far away from a similar solution, whereas in the quasiperiodic tiling, similar environments 
are guaranteed to occur within a much shorter distance - about twice the linear size of 
the region - thus allowing a greater overlap. The quasicrystal wavefunctions are probably 
"critical" states, with a power law long distance behavior and truly localized states should 
exist only exceptionally (such as the ring- states for special energies mentioned further 
below). Although one dimensional quasiperiodic systems have been successfully tackled by 
renormalization group methods Hamiltonians in two and higher dimensions have been 
mainly examined numerically, for example in the 2D Penrose and octagonal tilings 

It is known in the case of disordered Hamiltonians that level spacing statistics and the un- 
derlying wavefunctions fall into two categories : Wigner-Dyson gaussian statistics associated 
with extended electronic wavefunctions in the case of weak disorder, and Poissonian statis- 
tics associated with exponentially localized wavefunctions in the case of strong disorder. 
Although quasicrystals are not necessarily disordered, the complexity of the Hamiltonian 
makes it necessary and useful to carry out such statistical analyses, perhaps to discover a 
new type of statistical behavior in these media. 

In the work now to be discussed, we have used the inflation symmetry of quasiperi- 
odic systems to define a set of recursively related probability distributions for level statistics. 
Inflation refers to the (reversible) operation in which a subset of vertices of the tiling are 
erased, and the new vertices reconnected according to a well-defined prescription. The 
infinite quasicrystal is left unchanged by this transformation, aside from the trivial multi- 
plicative change of edge lengths by a factor A and a corresponding reduction in the vertex 
density by A*^ in d dimensions. For the octagonal tiling, this factor is A = 1 -|- \/2. This 
symmetry was used in in a discussion of the spectrum of tight-binding models on the oc- 
tagonal tiling, a full renormalization scheme (such as in for the Fibonacci chain) has not 
yet been obtained for the 2D tilings. Previous numerical studies did not take this property 
of inflation into consideration in analysing level statistics. However inflation symmetry is 
an important symmetry in the absence of translational invariance and the recursive scheme 
proposed in this paper is an attempt to utilise this property in solving for the level spacing 
probability distribution. 

The systems that we consider are square pieces of tiling whose interiors reproduce the 
quasiperiodic arrangement found in the perfect infinite octagonal tiling, and with edges that 
permit the periodic repetition along the two directions of the plane. These "square ap- 
proximants" of the octagonal tiling can be obtained by projection from a four dimensional 
cubic lattice by the "cut-and-project" method P]. Fig. la shows the k = 3 approximant, 
containing 239 sites, while larger square approximants are generated using the series of ra- 
tional approximants: A^^ = |,f,^---- Here we present results obtained for k = 4,5,6, 
containing 1393, 8119 and 47321 sites respectively (the number of sites A^^ increases by ap- 
proximately A^ with each increase in k). Successive approximants are related to each other 
by an inflation operation. The tight-binding Hamiltonian is defined by taking a hopping 
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matrix element t = 1 between connected sites (see Fig.l), and imposing periodic boundary 
conditions. The resulting sparse matrices were diagonalized by Lanczos routines to obtain 
the energy levels. Because of the presence of a discrete symmetry (reflection about the line 
x=y) it turns out that adding a Bohm-Aharonov flux does not change the model from the 
orthogonal to unitary symmetry class 0. This fact was used to augment the statistics for 
the smaller approximants by varying the boundary conditions and obtaining a large number 
of spectra. The level spacings are calculated independently for each subspace, we present 
here the results for the symmetric subspace (even parity under x- and y- translations and 
reflection across the diagonal), with 21750 nondegenerate levels for the biggest size consid- 
ered. Generally speaking, all the levels are in fact non-degenerate, aside from a macroscopic 
degeneracy of the E = states in this model, due to the eightfold ring-states that are 
permitted in this geometry - Fig.l illustrates such a wavefunction: localized on the eight 
highlighted sites whose wavefunctions are of constant magnitude but alternating sign - this 
degeneracy is well-understood and does not affect the level statistics. 

Fig. 2) shows a typical density of states (DOS) histogram obtained for the approximants. 
Its characteristically spiky shape poses the problem of "unfolding" the level spacings: in 
conventional random systems, fluctuations of the density of states (DOS) remain small 
enough with respect to a more smoothly varying underlying contribution, which is then 
used in renormalizing or unfolding the level spacings. In the randomized versions of this 
model that were considered in |^ this unfolding is possible as the density of states has 
an underlying smooth part. But in the absence of randomness, the DOS has no smooth 
variations, and the limit limN^ooP(E) for fixed E does not exist: the "fluctuations" are of 
the same order as the "average" DOS and they are present on arbitrarily small energy scales. 
We consider first therefore the "bare" (i.e. sans unfolding) level statistics for these systems. 
Fig. 3a shows the numerically calculated points and the fitted Lognormal (LN) distribution 
first reported in for the biggest approximant and Fig. 3b shows the points obtained for 
three different system sizes, after shifts so as to be centered at the origin. The points lie on 
gaussians of the form Pfc(lns) = ,^y^Q-°'i^^^-"-k) _ ^he peak position shifts with increasing 
k, reflecting the fact that the mean value of s is reduced by A^. More surprisingly, the width 
does not seem to depend very much on system size. The fitted values a vary from about 
0.79 to 0.81 from the smallest to largest system size, with an uncertainty of about 0.01 due 
to the fluctuations and the deviation from gaussian behavior in the tails. The continuous 
curve in Fig. 3b is obtained from the recursive calculation discussed further below. 

LN distributions may be obtained from iterative schemes such as the following: at each 
step of the iteration intervals are split into two nonequivalent subintervals with some choice 
of renormalisation factors. After n steps, for n sufficiently large, one approaches a LN dis- 
tribution of interval lengths as the saddlepoint approximation used becomes better (such 
a calculation was done in ^ for a resistive network). This heirarchical construction is ap- 
pealing for quasiperiodic systems, with their inflation-invariant structure and would predict 
that both the mean value as well as the variance of Ins are linear in n . However since this 
latter behavior is not observed , the distribution of spacings for the quasicrystal must be 
more compact than this simple model predicts, implying a more "rigid" spectrum. 

We now examine unfolding procedures to correct the level spacings for DOS variations. 
A first attempt was reported in |jlO[ where each spacing was divided by a local mean value 
of the level spacing. One can set s^-™'* = (-Ej+i — Ei) / { ^'^'^^'~"^ ) where the "most local" 
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unfolding corresponds to taking m = 1. The resulting variables {s^"^^^^} were found to follow 
a Wigner-Dyson distribution. As m is increased, the probability distribution of the variables 
shifts continuously over to the LN distribution. This result, that the level spacing obeys 
W-D statistics upon unfolding was rediscovered by Zhong et al |TT| where they smoothed 



the integrated DOS using a spline fit. This works because the IDOS, while not a smooth 
curve, has eigenvalues that are grouped together in a way as to allow a piecewise continuous 



fit - this procedure yields the same result as in |jTO[ for m small. 

To resume, the spacing statistics are W-D when the unfolding procedure is local, i.e. 
spacings are measured by the yardstick of the local DOS. At the other extreme, when the 
unfolding is done using the global mean level spacing the LN distribution is obtained. These 
unfolding schemes do not make use of the inflation transformation relating the spectra of 
tilings of successive sizes. With each inflation, there are times more levels eI''\ which fall 
within the energy intervals of the smaller system in an inhomogeneous fashion. Otherwise 
put, each "ancestor" interval Si of the {k — l)th approximant splits into a certain number 
rii of intervals of the (larger) kth approximant, with the mean value of Ui being equal to 
A^. We propose thus to renormalise each level spacing of the larger system by a mean local 
spacing, deflned with respect to the common ancestor as follows: 

(fc) p(fc) p(fe) 

where the indices i and j run from 1 to A^^ — 1 and A^^-i ~ 1- ^ indexes a given spacing 
of the bigger system, while j is the index of the ancestor spacing which has split into Uj 
subintervals after inflation, to the {k — l)th approximant. Fig. 4 shows the probability 
distribution of the new variables x-'^^ for the case k = 5 along with the Wigner-Dyson 
curve Pwd{x) = ^xe~'^^ 1'^. We note that the x variables as deflned automatically have 
the mean value of unity, and that no 'parameters are involved in the flt to W-D statistics. 
For the square approximants, one can thus deflne variables that correspond to the unfolded 
level spacings and which have Wigner-Dyson statistics independently of the system size (for 
systems large enough so as to provide reasonable statistics). The flt to Pwd is not as good 
as in [lTO| , pA| , due to the fact that the transformation A; — > A; + 1 is a discrete transformation 
- if it were possible to tune A gradually up from unity to its flnal value, one would have 
an unfolding that was truly local at each step and the resulting x— distributions would be 
exactly Pwd- 

Consider Eq.|] after one more iteration: 

(fe) 

(fc) (fc-l) Si /„x 

x\ 'x) = njm^j^ (2) 

where sf^^ is a subinterval of Sj^'^^ (subdivision index nj), which in turn is a subinterval of 

(subdivision index rii). The x variables on the kth and the {k — l)th levels obey the 
WD statistics. 

The quantities that contain information about this particular quasiperiodic system are 
the subdivision factors Uj. The distribution of Uj is not a smooth function as the variable 
takes on only a very limited set of values. However, when the inflation is carried out twice, 
the resulting subdivisions have a wider spread of values and a new distribution can in fact 
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be defined. The distribution of the product of subdivision factors appearing in Eq.^ was 
calculated for our set of three tilings, with the result shown in Fig. 5 {n is the product 
variable). We propose the functional form Qn^ = i3'^ne~^"', for the analytical computation 
that follows, where (3~^ = ~ 34 (shown as a continuous line in Fig. 5). This implies 
that the distribution of splitting factors follows a Poisson distribution when — > + 2 as 
if levels are randomly distributed within the bandwidth. Although a theoretical calculation 
for this distribution is lacking, one may speculate that it probably holds for other inflation 
symmetric systems as well, with an appropriate choice of the parameter f3. Note that the 
fluctuations around the smooth curve in Fig. 5 are still rather large - whereas for higher 
products n = U^^-^rii {p > 2) more continuous curves will be had, enabling better fits of the 
corresponding distributions Qn^ . 

According to Eq.0 the level spacings at the kth and {k — 2)th level approximants are 
related by Ins^'''^ — Ins^''"^^ = in(^^^^^^^^— ^) = Y. The distribution function of Y, -Py(l^), is 
thus given by the convolution of the distribution functions of Ins^'^^ and Ins^*^"^-* : in terms of 
the Fourier transforms, -Py(Q') = P^^\q) P'^''~'^\—q) . We find the distribution of the product 
^(fc)^(A;-i)^^ = X = e^ by first integrating over the x variables, finding 



X 



est / dn n'^ g(f) K^{-kX/2) (3) 



(2) 

where Kq is the zeroth order Bessel function. Putting in the explicit form of Qn , one 
obtains 

Py = cst -^-^ F(4, (4) 
(1+X)4 ^ ' 2' 2' l+X^ ^ ^ 

where F is the hypergeometric function {X = 7rX/(2/5)). This function is very well approx- 
imated by a gaussian form for Py{Y), for Y close to its most probable value Y^ax ~ —4.1. 
One finds upon expanding around Y^ax that 

PY{Y) = cst e-"(^-^— (5) 

where a = 0.42. This along with the convolution relation Py = p(^^2)p('^) leads us to 
propose a gaussian form for the level spacings, 

P^'^) (Ins) = est Q-Ml^s~{^r,s),f ^g^ 

with a similar relation with k replaced by /c — 2 for p(^~2). From the convolution relation 
for Py one has 

PY{Y) = cst e~ki»^+^^-2)iy-{y)r (7) 

where (Y) = (Ins)^ — (lns)/c_2- Comparing Eqs.^ and|^ for Py, one gets 

^{ak + ak-2) = o- ; (lns)fc_2 - (Ins)fc = y^a^ (8) 

Thus if a is independent of k, one has a = 0.84, the value taken for the gaussian curve 
plotted in Fig.3b. The numerical data shown in Figs. 3 are in accord with this, if one notes 
that the fitted gaussian width increases slightly with system size and approaches the value 
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given above. The shifts calculated numerically were (lns)fc+i — (Ins)^ = 1.8 for k = 4 and 
k = 5, while the expected value from the calculation above is ^InYmax ~ 2. 

The LN distribution is clearly only valid for values of Ins close enough to the mean value: 
the deviations from the gaussian approximation to the exact expression in Eq.|^ are most 
evident in the small Ins tail of Fig. 3a. 

In conclusion, level statistics in our quasiperiodic system are described by both Wigner- 
Dyson statistics and the LN statistics depending on the choice of unfolding used. The former 
indicates that the model belongs in the class of nonintegrable models described by random 
matrix theory, while the latter has its origin in the inflation-symmetry of the tiling. It would 
be interesting to consider related models such as the Penrose tiling (with inflation parameter 
r = ^~'^ ), where similar results should hold with numerical differences only in the width 
and the peak shifts of the LN distribution of spacings. 
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FIGURES 




FIG. 1. Square approximant containing 239 sites, showing the eight- fold symmetric sites 
around which the E = wavefunctions are locahsed (see text). 
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FIG. 2. Histogram of the density of states of a square approximant 
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FIG. 6. Plot of the distribution of values of the product TLin.j of subdivision indices in going 
from the {k — 2)th to kih approximant. The curve is the analytical form discussed in the text. 
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